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Abstract 

We determine the critical momenta for chiral fermions in the domain wall model re- 
cently suggested by Kaplan. For a wide range of domain wall masses m and Wilson 
couplings r we explicitly exhibit the regions in momentum space where the fermions 
are chiral. We compare the critical momenta for the infinitely large system with those 
obtained on a finite lattice. 



1 Introduction 



In a recent paper a new method for simulating chiral fermions on a lattice was suggested. 
The proposal uses the fact that by introducing a domain wall in an odd dimensional - 
and therefore vectorlike- theory one finds zero modes bound to this domain wall 0, ^, 
Regarding the domain wall as a lower (even) dimensional world it was demonstrated in []T[ 
that for the infinite lattice these zero modes are chiral fermions and that the resulting theory 
on the domain wall exhibits the desired anomaly structure. 

The same features have been shown to survive on a finite lattice where boundary condi- 
tions require one to consider a wall/anti-wall pair. For low energies a chiral fermion is bound 
to one of the domain walls with its chiral partner of opposite chirality living on the other 
domain wall. In addition, it has been demonstrated that in the presence of weak external 
gauge fields the divergence of the gauge current satisfies the continuum anomaly equation 
. This surprising result has been explained as being a consequence of the fact that on the 
lattice the divergence of the Goldstone-Wilczek current |^ has the same form as in the 
continuum 0. 

One of the most important properties of the domain wall method is that the fermions 
are only chiral in the low energy limit. For the special choice of the Wilson coupling r = 1 it 
was shown in that there exists a critical value of the momentum pc for which the fermions 
cease to be chiral. However, it is expected that also for other choices of the Wilson coupling 
such critical momenta exist. This was confirmed on the finite lattice where the existence of 
the critical momentum could be demonstrated for an r = 1.8 ||^. 

Here we want to perform a systematic study of the values of the critical momenta for 
which one will find chiral modes. We compare the analytical results obtained on the infinite 
lattice with the spectrum of the finite lattice Hamiltonian. Aside from the theoretical interest 
of this question, our results are of practical importance for future numerical simulations of 
this system. 



2 Zeromodes on the Infinite Lattice 

To be specific we will first discuss the zeromodes for the case of a 3-dimensional model, 
though our results will be generalized to arbitrary dimensions at the end. We start with the 
Dirac- Wilson operator on an infinite lattice with lattice spacing a = 1 

3 3 

KsD = E + me{s) + ^ E (1) 

where d denotes the lattice derivative d^ = \ \^z,z+ii ~ ^z,z-ij\ , ^ the lattice Laplacian = 
[^z,z+^l + Sz,z-fi — '^Sz,z] 5 the 0"^ are the usual Pauli matrices and r the Wilson coupling. We 
will denote by s the extra dimension along which the mass defect appears, while x, t are the 
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2-dimensional coordinates. The domain wall is taken to be a step function 6, 

-1 s < 

s = (2) 

+ 1 s > 

where the height of the domain wall is given by the mass parameter m which we will choose 
to be positive throughout this paper. 

We are looking for solutions which are plane waves in the (x,t)-plane 

^± = e^(pt*+P-^)$(s)^± (3) 

where are the eigenspinors of o"3 

a^iU^ = ±u^. (4) 
With this ansatz the Dirac operator becomes 

2 2 ^ 

KsD = sin(pi) + a^ds + m6'(s) + r(E(cos(pj) - 1) + - A, . (5) 

1=1 i=l 

Our final goal is to diagonalize the 3 dimensional Dirac operator in such a way that it reduces 
to the 2 dimensional Dirac operator for free massless fermions, 

7^3^^ = K^d"^ (6) 

where K2D acting on \1' is given by 

2 

K2D = '^f^^f^ = sin(pt) + (T2 sin(p^)) . (7) 
^=1 



Hence the equation to solve is 

' $M± = (8) 



asds + meis) -rF + ^A, 



where 

F= 5:(l-cosfe)). (9) 

i=t,x 

Following |1| we choose an exponential ansatz for $ away from the domain wall 

+ 1) = 2<l>(s). (10) 
Inserting this into (H) one finds four solutions 



r - rueff ± ^TUeff {nieff - 2r ) + 1 

z = — ^ 11 

r ± 1 
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where meff = m6{s) — rF, the ± in the nominator stand for the two roots and the it in the 
denominator stand for the chirahty. Note that in the hmit r = 1 eq.(|ll]) can be reduced to 
the corresponding expressions in [|I[|. 

We have to impose the condition that the solutions are normahzable to obtain sensible 
wavef unctions. This means that |z| > 1 for s < and |z| < 1 for s > 0. The boundaries of 
the regions where chiral solutions exist are obtained by setting \z\ = 1. Explicit matching 
of the normalizable solutions for positive and negative s at s = enables us to determine 
the regions with chiral fermions. We find that existence and chirality of the solutions is 
independent of the sign of r and that a negative m leads to opposite chiralities. Depending 
on the values of m/r we get m = rF and m = r{F + 2) as the boundaries for the critical 
momenta where F is defined as above in eq.@. 

The results are summarized in fig.l. We show the Brillouin zone — vr < < vr, 
— TT < < TT for different ratios of m/r. The white areas indicate the region in momen- 
tum space where chiral fermions exist. Starting with m/r = we find no chiral fermions. 
For increasing m/r > the region where chiral modes exist grows from a small circle around 
p = (0,0) until it hits the boundary of the Brillouin zone for m/r = 2. The boundary of 
the circle, i.e. the upper critical momenta, is given by m = rF. Increasing m/r further 
opens up the two "doubler" modes at p = (0, vr) and p = (vr, 0) which have fiipped chirality, 
while the original mode at p = (0, 0) disappears. Here the boundaries of the white regions 
are given by m = rF for the lower and m = r{F + 2) for the upper critical momenta. For 
m/r > 4 the two "doublers" disappear and we get a zero mode at p= (vr, vr) with the same 
chirality as the mode at p= (0, 0). The boundary for the lower critical momenta is given by 
m = r{F + 2). This mode is finally also lost as m/r is increased to m/r > 6. 

We want to remark that this spectrum stems from solutions only, and that there are 
no ^~ solutions for positive m. The change of the chirality is the usual reinterpretation of 
the chirality at different corners of the Brillouin zone. 

The generalization to arbitrary dimensions d = 2n + 1 consists merely in replacing the 
function F in eq.(^) by 

F = U(l - cos(p,)) . (12) 

i 

We find in d dimensions that for 2k < \m/r\ < 2k + 2 {0 < k < d — 1), the number of chiral 
zero modes Nzm bound to the d — 1 dimensional domain wall is given by 

iV.r.= ('^;^) (13) 

and their chirality is (—1)'^. 

The regions of chiral fermions as found here correspond exactly to the values of m/r 
where the lattice Chern-Simons current induced by Wilson fermions changes its value 
giving contributions to the Goldstone-Wilczek current. The calculation of the Chern-Simons 
current in uses the observation that the fermion propagator in momentum space can be 
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interpreted as a map from the torus to the sphere {T'^ —>■ S'^). The winding number of this 
map is closely related to the zeromode spectrum and changes only at particular values of m/r 
which agree exactly with our results for the values of m/r where the number of zeromodes 
changes. 



3 Finite Size Effects 



Any numerical work on this system will necessarily involve finite lattices, and so we now 
compare the critical momenta obtained on the infinite system with the ones of a finite 
lattice. On the finite lattice we have to choose some boundary conditions which generate a 
second anti-domain wall. The mass term is therefore modified to be itlOl^s) with 

r -1 2 < s < ^ 
Ol{s) = \ +1 ^f + 2<s<Ls . (14) 
[ s = 1,^ + 1 



We searched for the zeromodes on the finite lattice by solving the Hamiltonian problem 
numerically. If we again assume plane waves in the x-direction the Hamiltonian is given by 



H = -ai 



r 

sin(pa;) + (y'ids + m6L{s) + r{cos{p^) - 1) + -A^ 



(15) 



We computed the eigenvalues and eigenfunctions of the Hamiltonian eq.(|T3]) numeri- 
cally 0. Choosing the Hamiltonian instead of the Dirac operator reduces the numerical effort 
substantially as we only have to search for a single critical momentum instead of a pair. To 
find the critical momenta we studied the ratio 

which is a normalized measure for whether the fermions are chiral or not. It is zero if the 
fermions are chiral and R> for non-chiral modes (see fig. 2b in [§). To determine whether 
we still have chiral fermions we defined a threshold value for R. If < 0.01 we regarded the 
fermions to be chiral. 

We show the critical momenta as functions of r at three different values of m in Fig. 2. 
As discussed above we should find two boundary curves for the critical momenta. Note that 
because we now use the Hamiltonian, the function F is only one dimensional F = 1 — cos(pa;). 
For m/r < 2 the chiral fermions appear for momenta bounded by m = rF. This is the solid 

^In principle it is, of course, possible to perform a similar analysis as in the infinite system. However, 
on the finite lattice one also has to take into account solutions which are non-normalizable in the infinite 
system. In addition, on the finite lattice the modes are not exactly chiral. A mode bound to one of the 
walls always gets an exponentially suppressed contribution from the mode on the other wall. These features 
render analytic calculations on the finite lattice quite intractable. 



4 



line in Fig. 2. For 2 < m/r < 4 the momenta for which chiral fermions appear are given by 
m — rF for the lower and m — r{F + 2) for the upper critical momenta. We plot the curve 
for the upper critical momenta as a dashed line in Fig. 2. 

We compare the results from the infinite system with our finite lattice calculations. The 
crosses correspond to a system size of L = 100 and the open circles to a size of L = 20. We 
find that the L = 100 lattice is practically indistinguishable from the infinite system. For 
L = 20, a lattice size realistic for simulations, a small shift occurs. Fixing m and r we find for 
m/r < 2 a smaller value and for 2 < m/r < 4 a larger value of the critical momentum. Note, 
that for m = 0.4 the circles belong only to the solid curve. We did not find the momenta 
which correspond to the dashed line as our resolution in the numerical computation was not 
fine enough. 

In conclusion, we find the differences between L = oo and L = 20 to be small. This 
means that the lattice has not to be too large to reproduce the basic features of the model at 
L = oo which makes the domain wall model feasible for numerical investigations. Therefore 
we have shown that the domain wall model can be used for numerical simulations on realistic 
lattices. We also give the values of the domain wall mass m and the Wilson coupling r with 
which numerical simulations should eventually be performed. 
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Figure Caption 



Fig.l We plot the regions within the Brillouin zone — tt < pt < +7r, — tt < Px < +7r where 
chiral fermions exist (white areas) as a function oim/r with m the domain wall mass and r 
the Wilson coupling. For the different cases < m/r < 2, 2 < m/r < A and 4 < m/r < 6 
we have one chiral fermion with positive chirality, two chiral fermions with negative chirality 
and again one chiral fermion with postive chirality, respectively. For m/r < and m/r > 6 
there exist no chiral fermions. 

Fig. 2 We plot the two lines which give the upper and lower critical momenta on the infinite 
system, m = rF (solid lines) and m = r{F+2) dashed lines, where F = l~cos{px). We show 
these lines at three different values of the domain wall mass m as a function of the Wilson 
coupling r. We compare the curves from the infinite system with results from finite lattice 
calculations with lattice sizes L — 100 (crosses) and L — 20 (open circles). For m = 0.4 our 
resolution in the numerical computations was not fine enough to find the critical momenta 
corresponding to the dashed line. Note that though a shift in the critical momenta on the 
finite lattice is visible we find the same structure as for the infinite system. 
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